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Abstract. In this paper, we investigate some chemokinematical properties of the Milky Way disk, by using a sample 
composed by 424 late-type dwarfs. We show that the velocity dispersion of a stellar group correlates with the age of 
this group, according to a law proportional to i ' 26 , where t is the age of the stellar group. The temporal evolution 
of the vertex deviation is considered in detail. It is shown that the vertex deviation does not seem to depend 
strongly on the age of the stellar group. Previous studies in the literature seem to not have found it due to the use 
of statistical ages for stellar groups, rather than individual ages. The possibility to use the orbital parameters of 
a star to derive information about its birthplace is investigated, and we show that the mean galactocentric radius 
is likely to be the most reliable stellar birthplace indicator. However, this information cannot be presently used 
to derive radial evolutionary constraints, due to an intrinsic bias present in all samples constructed from nearby 
stars. An extensive discussion of the secular and stochastic heating mechanisms commonly invoked to explain the 
age-velocity dispersion relation is presented. We suggest that the age-velocity dispersion relation could reflect the 
gradual decrease in the turbulent velocity dispersion from which disk stars form, a suggestion originally made by 
Tinsley & Larson (19781 and supported by several more recent disk evolution calculations. A test to distinguish 
between the two types of models using high-redshift galaxies is proposed. 
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1. Introduction 

The publication of the data from HIPPARCOS (ESA 
I1997fl improved substantially kinematical studies of the 
solar neighbourhood. Simultaneous measurements of par- 
allaxes and proper motions, made by the astrometric satel- 
lite, allow the calculation of tangential velocities for tens 
of thousands of stars. In spite of it, a complete study of 
the stellar velocity distribution in the vicinity of the Sun 
is yet to be made, since the database for radial velocities 
has not increased in the same proportion, and is still a 
major problem in stellar kinematic studies. 

The correlation between kinematical and chemical 
properties of the disk stars is a classical result of galac- 
tic astronomy (Roman 119501 HUS^I Eggen, Lynden-Bell & 
Sandage 1962, Edvardsson et al. 119931 hereafter Edv93). 
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The knowledge about both the kinematical and chemi- 
cal properties of the stellar populations gives several clues 
about the evolution of a galaxy, unveiling the complex 
processes that concurred in its building. 

Recently, we have built a large sample having chromo- 
spheric ages and photometric metallicities. This sample 
was used to study the age-metallicity relation (Rocha- 
Pinto et al. BODOaJ hereafer Paper I) and the star forma- 
tion history (Rocha-Pinto et al. I2000bl hereafter Paper 
II). It is interesting to use these data in the study of 
chemodynamical constraints. Our sample is small, com- 
pared to that used in some recent investigations, but it 
can still be useful for the investigation of some impor- 
tant constraints. Several papers that address constraints 
between stellar ages and kinematical properties make use 
of ages calculated by statistical methods, that is, average 
ages for a group of stars. The advantage of our approach is 
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that we consider absolute ages, and the constraints are ex- 
pected to give a more accurate description of the galactic 
kinematics. 

This paper is organized as follows. The kinemati- 
cal sample is described in detail in Section 2. The age- 
velocity dispersion relation (hereafter, AVR) is derived in 
Section 3, where we show that all components of the spa- 
tial velocity of the star present a similar AVR. In Section 4, 
we investigate the temporal evolution of the vertex devi- 
ation. Section 5 investigates the possibility of using or- 
bital parameters to find the stellar birthplace and age. In 
Section 6, we attempt to use orbital information to derive 
some radial constraints to chemical evolution parameters. 
The AVR is readdressed in Section 7, where we present a 
thorough discussion of its interpretation, concluding that 
secular disk heating mechanisms are incapable of account- 
ing for the large dispersions of old disk stars, and suggest- 
ing that a major component of the AVR may be due to 
disk cooling, as originally suggested by Tinsley & Larson 
(1978). Our conclusions are presented in Section 8. 

2. The kinematic sample 

We have looked for kinematical data in the literature, 
for the 729 stars of the initial sample and for the 261 of 
the aditional sample (see Paper I for details). Proper mo- 
tions were obtained from the HIPPARCOS catalogue or 
the SIMBAD database. Radial velocities come from sev- 
eral sources: Barbier-Brossat, Petit & Figon (1994), Duflot 
et al. (|1995al I1995b|> . Fehrenbach et al. ifuMl HW^i . 
Nordstrom et al. H1997[> and in the SIMBAD database. 
In total, 459 stars have radial velocities in the literature, 
from which 424 have chromosphcric ages lower than 15 
Gyr. In the following, we will use only these 424 stars. 

We calculated the spatial velocities U, V and W, in 
a system pointing to the galactic anticenter, the direc- 
tion of the galactic rotation and the direction towards 
the north galactic pole, respectively, using the equations 
of Johnson & Soderblom (1987). The velocities are fur- 
ther transformed into peculiar velocities with respect to 
the local standard of rest, u, v and w, using the basic 
solar motion 1 (Mihalas & Binney 1981), (uq,Vq,Wq) = 
(-9,ll,6)kms- 1 . 

The sample is given in TableQ] available at the CDS 2 , 
where we list, respectively, the HD number (or BD num- 
ber for a few stars), the chromospheric index logi?^ K , 
the photometric metallicity [Fe/H], the radial velocity, the 
chromospheric age, calculated according to the formalism 
described in Paper I, the peculiar velocities u, v and w 
with their measurement errors, and some parameters re- 
ferred to later on. 

Spatial velocity distributions are strongly dependent 
on the mean age of the stellar group studied. We can con- 

1 The value of vq is still uncertain. Dehnen & Binney ( 1998 ) 
have calculated vq — 5.25 km/s, but Feast (2000) argues that 
several recent determinations give vq ~ 11 km/s. We decided 
to use this more conservative value in our analysis. 
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sider stellar subgroups having different ages, to investi- 
gate their velocity distribution. The most direct division 
in our sample is set by the chromospheric activity lev- 
els, since it was constructed from chromospheric activity 
surveys. Active and inactive stars must present different 
velocity distributions, that will reflect the mean ages of 
these subgroups. The active stars are expected to be more 
concentrated in the velocity space, compared to the in- 
active stars, as shown by Soderblom (1990) for a sample 
of 32 active stars. The same result was found by Jeffries 
& Jewell (1993), from a sample of dwarfs detected in the 
extreme ultraviolet by ROSAT. 

In Fig.^ we show the velocity diagrams for these two 
subgroups of chromospheric activity. The diagrams for the 
active stars are presented in the left panels, and those for 
the inactive stars in the right panels. The plot confirms the 
previous conclusions. Active stars present a low velocity 
dispersion, that is explained by their youth, compared to 
the inactive stars. Note that some active stars present ve- 
locity components typical of an advanced age. Soderblom 
( 199(3) nas already pointed to this problem. However, in 
his sample of 32 dwarfs, only 3 showed high velocity com- 
ponents. Rocha-Pinto & Maciel l(T§9"gl hereafter RPM98) 
have increased this number to 10. The present sample 
contains 29 stars like these. Jahreifi, Fuchs & Wielen 
(1999) also mention the existence of chromospherically 
active high velocity stars. Comparing isochrone, chromo- 
spheric and lithium ages of these objects, Rocha-Pinto et 
al. ( 2002) concluded that they can be single stars formed 
out of the coalescence of short-period binaries. The coa- 
lescence makes the star look chromospherically younger, 
because the orbital angular momentum of the former bi- 
nary is transformed into rotational angular momentum of 
the newly formed single star; however, since the binary 
formed some Gyr before the coalescence, the movement of 
its center of mass around the Galactic center, later inher- 
ited by the single star, shows the kinematical signatures 
of an advanced age. We call these stars CYKOS (chromo- 
spherically young, kinematically old stars). 

3. Age-velocity dispersion relation 

The velocity dispersion of a stellar group is larger, the 
older is the group. This relation is usually interpreted as 
a result of gravitational perturbations experienced by the 
stars during their translations around the galactic center. 
In spite of the fact that several processes which can cause 
an increase in the velocity dispersion are relatively well 
known, it is not completely understood which of them is 
the main responsible mechanism in this heating, and this 
comes mainly from the disagreement between the observed 
dispersions amongst disk dwarfs found by different stud- 
ies. Several authors have studied this issue fWielen 119741 
Mavor 11371 Cayrel de Strobel 11371 Carlberg et al. 11381)1 
Knude et al. 111)871 Stromgren lT3571 Meusinger et al. lTMTl 
Wielen et al. H"§9"2l Dehnen & Binney IH)9"8l Fuchs et al. 
2001) but their results are quite conflicting. It is not clear 
whether this was caused by selection effects, a problem 
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Fig. 2. Spatial velocities as a function of age. The disk 
heating can be clearly seen in this picture. Some stars, 
marked in the plot by solid circles and empty triangles, 
were discarded from the sample, before the final calcula- 
tions of the age-velocity dispersion relation (see text). 



which plagues severely kinematical studies, or is an in- 
trinsic property of different stellar classes used by them 
(see Fridman, Khoruzhii & Piskunov ll994f) . 

In Fig. 12 we present plots of peculiar velocities as a 
function of the stellar age. It can be seen clearly that the 
number of high- velocity stars grows towards older ages. 
In the middle panel, v clearly shows the asymmetric drift 
(Binney & Merrifield 1998), that is, the older stars ro- 
tate around the galactic center with slower velocities with 
respect to the local standard of rest. 



In order to calculate the velocity dispersion we had to 
further examine our stellar sample. First of all, there are 
three clear outliers, as seen in Fig. [21 (solid circles). For 
these outliers at least one of the velocity components is 
not typical of thin disk stars. These are all CYKOS, which 
are not supposed to be real young stars, thus, we have 
discarded them from this analysis. Furthermore, we have 
also excluded from the velocity dispersion calculations all 
the stars with metallicity [Fe/H] < —0.5, since they are 
likely to belong to the thick disk. This condition excludes 
24 stars (empty triangles in Fig. 2), so we are left with 397 
stars (crosses) from the 424 star sample. 

Each component of the velocity dispersion is calculated 
in equal star number bins instead of equal age range bins. 
Because the sample stars are not evenly distributed in age 
we have used bins with equal number of stars (N stars / hin = 
40) . This way we have an equal weight for each data point 
when fitting the AVR. 

The AVRs that we have found can be approximated by 
power laws of the kind a oc (1 + t/r) x , where r is a time 
scale for the growth of the velocity dispersion and t is the 
stellar age. The exponent x is 0.26±0.01, 0.24±0.02, 0.30± 
0.03 and 0.34 ±0.02, for ctot, cr u , a v and a w , respectively. 
These exponents can be considered equal within 2a. The 
relation in w is slightly steeper. 

Nevertheless, this fitting is very sensitive to the last 
age bin. If it is removed, we find x ~ 0.31 for dtot- If, 
on the other hand, we use bins equally spaced in age, we 
find x equal to 042 ± 0.04, 0.40 ± 0.02, 0.48 ± 0.04 and 
0.41±0.03, for er tot , a u , a v and a w , respectively. The large 
difference between these exponents and those calculated 
for equal star number bins reinforces our concerns about 
the sensitivy of the observed AVR. The equal star number 
binning was chosen to avoid many of the uncertainties that 
arise from the use of age range bins with poor statistics. 

Other fitting laws commonly used in the literature 
were tried. For a parametrization like a = ao + a\t x , 
we have found x to t = 0.24 ± 0.04, x u = 0.23 ± 0.04, 
x v = 0.37 ± 0.07 and x w = 0.30 ± 0.09. For a oc t x , we 
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Fig. 1. Velocity diagrams v x u and v x w for the active stars (left panels) and inactive stars (right panels). Compare 
the increase in the velocity dispersion of the inactive stars with that of the active stars. Note that, in spite of their 
low velocity dispersion, amongst active stars there are objects with velocity components typical of an old population 
(Rocha-Pinto et al.EDEJ- 



have x to t = 0.26 ± 0.03, x u = 0.23 ± 0.04, x v = 0.28 ± 0.04 
and x w = 0.35 ± 0.06. 

Our age-velocity dispersion relations are given in Fig. 
13 which shows that our dispersions are considerably 
higher than those found by Knude et al. ( 1987), Stromgren 
(1987) and Meusinger et al. ((TMTj) . for the older ages. The 
agreement is substantially better for the earlier studies by 
Wielen fT374) and Cayrel de Strobel fl$7fy . The data 
given by Fuchs et al. (2001 ) follow closely those of Wielen 
l|1974[) . and are not shown in the figure. 

The data in Fig.Olare presented also in TableEl In this 
table, the first column gives the age range considered, and 
the other columns give, respectively, Otot and the error 
in this quantity; <j u , a v , a w , and their respective errors, 
each shown after the corresponding velocity dispersion; 
the axial ratios a v : a u and a w : a u ; and the asymmetric 
drift according to the age of the stellar group. 

Although our data agrees at first glance with those 
given by Fuchs et al. (2001), the fitted velocity dispersion 
laws in their work and ours differ considerably. It is impor- 
tant to understand the cause of this difference. Fuchs et 
al. have used a sample of late-type stars based on chromo- 



spheric ages, like ours. But they have not taken in consid- 
eration the metallicity dependence of the chromospheric 
indices (Rocha-Pinto & Maciel I1998J) . when calculating 
the stellar ages. Moreover, they have weighted each star 
by the modulus of the vertical velocity, \W\, before cal- 
culating the velocity dispersions. The \W\ weighting was 
introduced by Wielen (|1974f) to take into account the prob- 
ability for the detection of each star in the local sample, 
given their vertical oscillation with respect to the Galactic 
plane. 

By disregarding the metallicity dependence of the 
chromospheric ages, Fuchs et al. f2QflI} are likely to mix 
hotter older stars with somewhat 'warm' intermediate-age 
stars. This could boost the velocity dispersion to high val- 
ues at 4-6 Gyr, resulting in a steeper AVR. This is illus- 
trated more properly in Fig. 0] for u. We have investi- 
gated this effect with our sample. The chromospheric age 
was recalculated using the calibration given by Donahue 
(1993), that is, without any correction for the metallicity 
of the star. The stars were binned by groups of 40, as in 
the former cases. For a oc t x , we have xtot = 0.33 ± 0.05, 
x u = 0.29 ± 0.05, x v = 0.39 ± 0.05 and x w = 0.37 ± 0.08. 
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Fig. 3. Characterization of the kinematical sample with a velocity ellipsoid. The panels a, b, c and d stand for the 
relations between the age and the velocity dispersions in the components u, v, w, and for the total space velocity 
(af ot = cr„ + cr^ + cr^), respectively. The velocity dispersions were calculated with respect to the principal axis of the 
velocity ellipsoid for each age group. Panels e and f show the evolution of the axial ratios a v /a u and a w /a u of the 
ellipsoid. Panel g shows the asymmetrical drift as a function of age. In these plots the symbols indicate: error bars, 
our data, calculated by binning the stars in intervals of equal star number (N = 40 stars), without applying the \W\ 
weights; semifilled triangles, Cayrel de Strobel 11974.1 letters W, Wielen 1974; letters K, Knude et al. 119871 letters S, 
Stromgren 1987; empty circles, Meusinger et al. 1991. 



This result shows that by disregarding the [Fe/H] de- 
pendence in the chromospheric ages, the exponents of 
the AVR have increased by 6% to 39% from the sup- 
posedly real values. This is an important result, since it 
can shed some light into the different slopes found by dif- 
ferent groups, using different star groups. As Fridman et 
al. I)1994|l point out, there are two main points of view 
in the literature for the AVR exponent: Wielen and col- 
leagues have found higher values than Knude et al. I|1987|l . 
Stromgren (1987) and Meusinger et al. ( 1991 ). The results 
by the Heidelberg group are largely based on chromo- 
spheric age estimates that disregard the metallicity depen- 
dence of the chromospheric ages, while the later studies are 
based on isochrone ages which are free of this systematic 
error. Since our results use a corrected cromospheric age 
scale, we should expect to find exponents closer to those 
found from isochrone age studies. 



The \W\ weighting also affects considerably the re- 
sulting AVR. For (i a f, and using the uncorrected 
chromospheric ages, as above, as in Fuchs et al. ((7001), 
we have found x tot = 0.50 ± 0.08, x u = 0.36 ± 0.09, 
x v = 0.49 ± 0.13 and x w = 0.53 ± 0.07. Qualitatively 
similar values are found when we use the \W\ weight- 
ing with the corrected ages and our preferred heating law 
a cx (1 + t/r) x : x to t = 0.50 ± 0.02, x u = 0.47 ± 0.02, 
x v = 0.57 ±0.04 and x w = 0.46 ±0.03. However, these fits 
are somewhat worse than the previous ones for which no 
vertical weight was applied. 

We have explored the extent to which subgiant and 
thick disk star contamination can affect the AVR. Because 
these stars are expected to come from a hotter Galactic 
component, their presence would surely inflate the velocity 
dispersions of the oldest bins, making the AVR artificially 
steeper. Although our adopted metallicity cut [Fe/H] < 
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Table 2. Velocity dispersions, axial ratios and asymmetrical drift. 



(Age(Gyr)) 


(Ttot 


(T u 




(J w 




a ■ a 


(V) 

\ * I 


0.22 


29 ±3 


25 ±4 


12 ±2 


9± 1 


0.47 ±0.08 


0.36 ±0.06 


-2 


0.55 


33 ±3 


28 ±3 


16 ±2 


7± 1 


0.57 ±0.06 


0.26 ±0.04 


-4 


1.17 


33 ±3 


23 ±2 


19 ±5 


13 ±2 


0.83 ±0.14 


0.55 ±0.09 


-2 


2.11 


50 ±4 


41 ±4 


22 ±2 


18 ±2 


0.54 ±0.09 


0.43 ± 0.07 


-16 


2.84 


57 ±6 


47 ±6 


24 ±3 


22 ±3 


0.52 ±0.06 


0.46 ± 0.08 


-26 


3.57 


52 ±3 


41 ±3 


25 ±4 


20 ±3 


0.63 ±0.09 


0.49 ±0.10 


-11 


4.51 


56 ±3 


48 ±4 


21 ±2 


18 ±2 


0.43 ±0.07 


0.39 ± 0.05 


-15 


5.90 


68 ±4 


52 ±6 


38 ±6 


23 ±3 


0.73 ± 0.09 


0.44 ± 0.08 


-13 


7.71 


65 ±4 


53 ±5 


30 ±3 


24 ±3 


0.57 ±0.09 


0.45 ± 0.07 


-17 


10.67 


69 ±4 


51 ±5 


36 ±4 


30 ±3 


0.69 ±0.10 


0.58 ± 0.09 


-38 



Table 3. Properties of the best-partition groups. 
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Fig. 4. Increase of the u peculiar velocity with age, for un- 
corrected (panel a) and corrected (panel b) chromospheric 
ages. The effect of neglecting the [Fc/H] dependence of the 
chromospheric age calibration is seen in panel a, where 
the mixing of old hotter stars with intermediate-age stars 
boosts the velocity dispersion around 4-6 Gyr, yielding 
a steeper AVR, which otherwise would present a slower 
increase (panel b). 

—0.5 already prevents our AVRs from being severely con- 
taminated by thick disk stars, a proper separation between 
the Galactic components can only be made with both 
kinematical and chemical data (Nemec & Nemec 
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Ojha et al. HMBI Chiba & Beers I2UUU|> . This is a classi- 
cal problem of mixture decomposition that can be reason- 
ably addressed with a multivariate finite mixture analy- 
sis. We have used the EMMIX program 3 (McLachlan et al. 
1999), which fits normal mixture models to multivariate 
data, using maximum likelihood through the Expectation- 
Maximization algorithm (Dempster et al. 1977). An ini- 
tial group configuration must be provided to the program, 
from which it looks for the best partition amongst the 
data. The best partition is not very much sensitive to this 
initial guessing, which can be either a random grouping 
or a pre-grouping using a K-means clustering (Hardigan 
1975|. The number of groups, n, among which the data 
will be partitioned must also be provided as input. Non- 
parametric tests (e.g., Tarter & Lock 1993) or physical 
considerations can be used to select n. Since we are mostly 
interested in singling out eventual thick-disk contamina- 
tors, we will assume n = 2. We applied the partitioning 
algorithm to our initial sample of 424 stars, in the four- 
dimensional space defined by [Fe/H], u, v and w. The av- 
erage properties of the two best-partitioned groups are 
shown in Table |21 The statistical decomposition yields 
a young, kinematically-cold, metal-rich group and a pre- 
dominantly old, kinematically-hot, metal-poor group. It 
is reasonable to identify this later group with a different, 
contaminant population, although its [Fe/H] is somewhat 
lower than typical thick disk metallicities, suggesting that 
old thin disk stars may have interloped this group. EMMIX 



The program EMMIX is freely distributed for non-commercial 
use at 

http : //www . maths .uq.edu.au/" gjm/emmix/emmix .html. 
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outputs the a posteriori probability, pi, of each star be- 
ing member of group i. To avoid thin-disk interlopers, 
we have considered a likely contaminator every star with 
Pu > 0.75. After discarding these stars from the sample, 
the AVR was calculated with and without the \ W\ weight- 
ing, using the law a oc t x . The resulting exponents can be 
found in Tabled together with a summary of all other fit- 
ted AVRs. We have found that the presence of thick-disk 
and subgiant contaminators in the sample can increase 
the exponent of the power-law AVR by up to 40%, in all 
velocity components. This caveat must be taken into ac- 
count in every AVR investigation, especially in face of the 
growing evidence supporting the singularity of the thick 
disk (Robin et al. 115551 Gratton et al. l2UUU|) . 

In Figs.|2l3 and f, we present the value of the ratio be- 
tween the semiaxes of the velocity ellipsoid, a v : <j u and 
a w : <r u . The values we have found are closer to those in 
Wielen (|1974[) and in Cayrel de Strobel l|1974|) . as well as 
those found by Dehnen l(T55g)l . for the 14369 HIPPARCOS 
stars. The axial ratios seem to be nearly constant, in op- 
position to the findings by Knude et al. ; 19K7li. according 
to which the young groups have a velocity ellipsoid with a 
substantially high a v : <r u , that decreases for old groups. 
This behaviour, as well as the very high values for a w : a u 
found by Meusinger et al. I|199 1 ) could reflect basically 
selection effects in the data samples. Finally, in panel g 
we show the asymmetric drift, which is the trend of the 
rotational velocity of a stellar population lagging behind 
more and more slowly the LSR, the larger its velocity dis- 
persions (Gomez & Mennessier 119771 Binney & Merrifield 

In spite of the apparent constancy of the velocity axial 
ratios, a sudden jump is seen in both a v : <r u and <r w : <r u 
at 1.2 Gyr ago, amongst our data as well as in the analysis 
of Wielen Q1974[l . Knude et al. Q1987[l and Meusinger et al. 
(1991). We presently do not have an explanation for this. 
If this feature is real, as these different analyses suggest, 
this could be an important constraint to the AVR theory. 

4. Vertex deviation 

The stars show a velocity distribution that can be char- 
acterized by a velocity ellipsoid (Schwarzschild 1907; 
Mihalas & Binney 1981), with semiaxes cr u i, a v \ and u w i, 
that correspond to the principal velocity dispersions of the 
stellar group. The semiaxes of the ellipsoid generally do 
not correspond to the coordinate axes used for the cal- 
culation of the spatial velocities (i.e., u, v and w), and 
some calculations must be done for the determination of 
these axes. In practice, it is known that the semiaxes a u \ 
and <T V \ are in the galactic plane (for instance, see Dehnen 
& Binney H998fl . and the problem is restricted to the de- 
termination of the angle ipv that indicates the galactic 
longitude towards which the principal major axis of the 
velocity ellipsoid points. This angle is called vertex devi- 
ation. Once determined, the peculiar velocity dispersions 
can be calculated with respect to the semiaxes of the ve- 
locity ellipsoid. 
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Fig. 5. Projection of the velocity ellipsoid of the stars 
younger than 1 Gyr onto the lif-plane. The tilted ellipsis 
with respect to the coordinate axes illustrates the vertex 
deviation of the data points. 



The vertex deviation tpy is defined by 



Vv 



l 



arctan ■ 



2al 



(1) 



where a\ v is the covariance between u and v. We have 
used a larger bin (each having 80 stars) to increase the 
significance of the results. 

Figure[S]presents the data points used for the youngest 
age bin used in this analysis. The dashed ellipsis is the 
projection of the velocity ellipsoid onto the uw-plane. The 
vertex deviation amongst these youngs stars is clearly seen 
by the tilt of the ellipsis with respect to the coordinate 
axes. 

The classical behaviour for the vertex deviation can 
be found in Delhaye l|19t>5[) . The vertex deviation is nega- 
tive (—50°) amongst O and B stars, but positive (~ 30°) 
amongst young dwarfs (mainly A type stars). Later dwarfs 
also present positive vertex deviations, but considerably 
smaller. The vertex deviation for very red and old stars 
is quoted as nearly zero. This relation between vertex de- 
viation and spectral type have led many investigations to 
assume an implicit relation between the age of the stellar 
group and its vertex deviation (e.g., Oblak^983 Oblak & 
Creze H^55jl . since later dwarfs have a higher average age 
than early dwarfs. 

Surprisingly, we do not find a similar behaviour in our 
data, as can be seen from Table Either using all the 



8 H.J. Rocha-Pinto et al.: Chemical enrichment and star formation III. 

Table 4. Fitted age-velocity dispersion relations. 



parametrization 


binned 
by 


age 
calibration 


\W\ 
weight 






x v 


Xw 


note 


a oc (1 + t/r) x 


number 


RPM98 


no 


0.27 ±0.01 


0.24 ± 0.02 


0.30 ±0.03 


0.34 ± 0.02 




a oc (1 + t/r) x 


number 


RPM98 


yes 


0.50 ±0.02 


0.47 ±0.02 


0.57 ±0.04 


0.46 ± 0.03 




(7 OC t 


age 


KrMyo 


no 


U.4z ± U.U4 


n /in x n no 
U.4U ± U.Uz 


U.4o ± U.U4 


n A i _i_ n oq 
U.41 ± U.Uo 




a xt x 


number 


RPM98 


no 


0.26 ±0.03 


0.23 ± 0.04 


0.28 ±0.04 


0.35 ± 0.06 




a xt x 


number 


Donahue (1998) 


no 


0.33 ±0.05 


0.29 ± 0.05 


0.39 ±0.05 


0.37 ±0.08 




a <xt x 


number 


Donahue (1998) 


yes 


0.50 ±0.08 


0.36 ±0.09 


0.49 ±0.13 


0.53 ±0.07 




a = an + a,\t x 


number 


RPM98 


no 


0.24 ±0.04 


0.23 ± 0.04 


0.37 ±0.07 


0.30 ± 0.09 




a oc t x 


age 


RPM98 


yes 


0.24 ±0.06 


0.24 ±0.07 


0.19 ±0.05 


0.25 ±0.03 


pn < 0.75 


a xt x 


age 


RPM98 


no 


0.24 ± 0.04 


0.25 ± 0.06 


0.28 ±0.04 


0.25 ± 0.03 


pn < 0.75 



stars in our sample or just those with vertical velocity 
smaller than 15 km/s, there is weak indication of an evo- 
lutive trend in the vertex deviation, if any. Moreover, the 
greatest vertex deviation is not seen amongst the youngest 
stars. 

Almost all of the former studies in the literature that 
tried to measure the dependence of the vertex deviation on 
the age used groups of stars with several ages. An average 
age for the group was calculated from the observed prop- 
erties of stars, and this average age was used in the deriva- 
tion of the function tpy{t). The main difference between 
this approach and ours is that the former can only mea- 
sure integrated kinematical properties of a stellar group 
(as a function of spectral type or colour). They are cumu- 
lative kinematical properties. In our approach, since we 
have ages for individual stars, we can consider the kine- 
matical properties in each age bin, in a differential ap- 
proach. 

There are few studies in the literature that do use a 
differential approach for the study of the vertex deviation. 
Their results give some support to our findings, regardless 
of their authors conclusion on the contrary. Byl (^974) has 
used evolutionary isochrones to calculate average stellar 
ages of nine stellar groups ranging from 0.38 to 10.50 Gyr 
old. The evolution of the vertex deviation amongst these 
groups does not show the implicit smooth trend expected 
from the classical theory, although much of the scatter 
could be blamed on the small number of stars in each of 
the nine groups. Wielen l|1974[l has also found high values 
for tpv at advanced ages, and a variation uncorrelated with 
age. Moreover, Wielen has also argued that the near zero 
vertex deviation found amongst late-type dwarfs could be 
a selection effect: Classic vertex deviation studies usually 
limited the samples to those stars with orbital eccentricity 
smaller than a given value, to avoid the contamination of 
the sample with high-velocity stars. By doing this, they 
sample the phase space within an ellipsis with semi-axes 
pointing to the main Galactic plane axes, thus biasing the 
sample. For old, late-type stars, which have large velocity 
dispersions, this bias would favour low ipv measurements. 
By using a less biasing selection criteria, \Jv? + v 2 < 40 



Table 5. Vertex deviation. 



(Age (Gyr)) Vv 



all stars 




0.39 ± 0.20 


14.8° 


1.64 ± 0.52 


16.0° 


3.20 ± 0.45 


20.0° 


5.20 ± 0.81 


11.3° 


9.16 ± 1.81 


10.8° 


stars with \ W\ < 15 km/s 




0.40 ± 0.20 


17.5° 


1.58 ± 0.51 


9.1° 


3.28 ± 0.44 


18.6° 


5.22 ± 0.88 


14.5° 


8.47 ± 1.34 


0.6° 


stars with e < 0.15 and \ W\ < 


15 km/s 


0.46 ± 0.18 


23.4° 


1.57 ± 0.51 


14.0° 


3.20 ± 0.45 


4.9° 


5.18 ± 0.76 


-13.3° 


9.23 ± 2.06 


9.5° 


stars with \/u 2 + v 2 < 40 km/s 


i and \W\ < 15 km/s 


0.39 ± 0.20 


16.4° 


1.45 ± 0.49 


19.6° 


3.26 ± 0.45 


-1.5° 


5.12 ± 1.68 


-9.2° 



km/s and IW^I < 15 km/s, Wielen showed that ipy m the 
range 10° to 30° could be found amongst stellar groups 
thought to have nearly zero deviation. We have tested 
both selection criteria with our sample. The vertex devia- 
tion amongst our stars using the classical selection criteria 
(orbital eccentricity e < 0.15 and \W\ < 15 km/s) and af- 
ter applying Wielen's selection criteria are given in Table 
These data were calculated within each of the original 
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bins, using only the stars that passed the selection crite- 
ria. Since the number of stars passing the Wielen criteria 
at the oldest age bin was very small, we combined the two 
oldest age bins. The result is somewhat close to the ex- 
pected bias pointed by Wielen. However, it is clear that 
the value ipy varies substantially depending on the selec- 
tion criteria used. Same order variations in the measured 
ipy have also been reported when the sample is limited by 
a coordinate range (Oblak 1983 Nunez & Figueras 1984). 

This variation could be caused by statistical errors in 
the calculation of ipy. If we consider only the error prop- 
agation from the velocity errors (taking as ~ ±2 km/s 
for each star), through a Monte Carlo simulation, we have 
e^v ~ 0.4°. But the error also depends on the Poisson 
noise due to sample size. We have run a bootstrap rou- 
tine to randomly select stars in each bin (so that some 
stars were counted twice or more and some were neglected 
from a given bin) and calculated the vertex deviation. 
This bootstrap selection was performed 1000 times. We 
have found that the average error in ipy is ^ 7.5°. This 
indicates that statistical errors alone cannot explain the 
variations we have found. 

Even non-differential analyses have not made a strong 
case for a smooth vertex deviation-age relation. The data 
compiled by Mayor : 1972; shows a large scatter in the 
value of ipy for each, adopted coeval, stellar groups. 
Particularly interesting is the fact that active HK dwarfs 
in his sample present a smaller vertex deviation (12°) than 
inactive HK dwarfs (20°), a result that bears resemblance 
with the entries of Table [5J Mayor concludes that large 
vertex deviations are observed amongst stars of all ages, 
provided that their orbital eccentricity is low. Recently, 
Ratnatunga & Upgren l)1997f) used a maximum-likelihood 
analysis in a sample of more than 700 Vyssotsky K and 
M dwarfs and showed that their vertex deviation is signif- 
icantly non-zero, amounting to 12° ± 4°. 

The existence of a smooth decreasing vertex deviation- 
age relation could provide an important clue about the 
cause of the vertex deviation. Theoretical work on this 
subject are generally split between those theories that pro- 
pose the vertex deviation is a reflection of the initial con- 
ditions, which we will call congenital scenario, and those, 
which explain it as the effect of perturbations in the local 
Galactic potential, here called perturbation scenario. 

Congenital scenarios proposed include the natural evo- 
lution of the original stellar configuration at birthplace 
taken as a spiral arm (Wooley 119701 Yuan 11971(1 . the 
launching of stars recently formed at the spiral arms with 
postshock velocity (Hilton & Bash 1982) or emerging from 
an expanding shell (Moreno, Alfaro & Franco 11999(1 . a 
temporary arrangement in the phase space due to the 
epicycles (Byl & Ovenden 119731 Bassino, Dessaunet & 
Muzzio 1986) or oscillations in the stellar orbits (House & 
Innanen 1975) and the superposition of macroscopic mo- 
tions over the Galactic rotation (Kato 1970). Perturbation 
scenarios address the effects of non-axisymmetry in the 
Galactic potential (Sanz & Catala [13571 Cubarsf [MJUJl 



and gravitational perturbations caused by the spiral arms 
(Mayor UnTJIinTJ. 

So far, none of these hypotheses can explain satis- 
factorily the vertex deviation phenomenon in all spectral 
classes. Some of the explanations are particularly suitable 
for the earlier stars, like the launching of recently formed 
stars by spiral arms or expanding shells, while other are 
more appropriate to explain the deviation found amongst 
late-type stars. 

A third class of explanation has gained momentum in 
the last years. The substantial improvement in the obser- 
vational techniques has allowed the investigation of the 
velocity phase space of the solar neighbourhood with un- 
precedented resolution. These studies have showed very 
convincingly that the velocity phase space is not homo- 
geneous, but punctuated by small scale structures and 
moving groups (e.g., Dehnen 1998; Skuljan, Hearnshaw 
& CottrellUHSS- Dehnen ifuMjl has shown that these 
moving groups themselves are the main responsible for 
the vertex deviation, and it is easy to understand this: 
Clumpy structures in the phase space puts a large weight 
on , yielding a larger ipy . This idea was first raised by 
Palous (1986 ), who argues that the Hyades and Sirius su- 
perclusters contaminate local samples, imprinting on their 
velocity distribution the vertex deviation. Moreno et al. 
(1999) present a similar reasoning for the explanation of 
the negative vertex deviation amongst OB stars. 

This could explain why the vertex deviation varies so 
widely depending on the selection criteria used and the 
coordinate range of the stars used in its determination. 
These different criteria would discard totally or partially 
some moving groups from the phase space, and this could 
induce to uncorrelated variations of <j\ v . In this case, the 
vertex deviation problem would be reduced to explaining 
the surviving of old moving groups and their peculiar ve- 
locities. A very original and promising idea was introduced 
by Dehnen 1(2000(1 and further developed by Fux 1(2001(1 
and Muhlbauer & Dehnen ( 2003J), according to which sev- 
eral features of the velocity phase space can be explained 
by the kinematic interaction between the outer disk and 
the Galactic central bar. 

5. Orbital parameters and stellar birthplaces 

From its spatial velocities and distance, it is possible to 
calculate the orbit of a star around the galactic center. 
The orbit is characterized by some parameters: perigalac- 
tic and apogalactic distance, R p and R a , average galacto- 
centric distance, R m , orbital eccentricity, e, and maximum 
height above the galactic plane, Z max . When a star is born 
from a molecular cloud, e and Z max are nearly zero. With 
time, due to the orbital diffusion, both parameters will 
grow, and the stellar orbit will become less circular. The 
orbits were integrated in a galactic potential, consisting 
of an exponential thin disk, a spherical bulge and a dark 
halo, as described in Flynn et al. 1(1996(1 . The value for the 
orbital parameters were calculated for 380 of the 459 stars 
in the kinematical sample, which correspond to the num- 
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Fig. 6. Comparison between two proposed birthplace in- 
dicators. The solid line shows the prediction by Wielen et 
al. (dSJ). 

ber of stars that follow the some additional criteria used 
in Table 2 of Paper I, that is, stars should be closer than 
80 pc to avoid substantial reddening in the photometric 
indices, brighter than 8.3 mag in V magnitude, less active 
than logi?^ K > —4.20 and have nominal chromospheric 
age lower than 15 Gyr. 

Let us consider the problem of the determination of the 
stellar birthplace. Wielen Ijl977|l claims that it is impos- 
sible to know with accuracy the galactocentric radius of 
the stellar birthplace from its orbit, due to the diffusion of 
stellar orbits in phase space caused by local fluctuations of 
the Galactic gravitational field. In each epicycle, the star 
would be scattered in random directions, and this would 
destroy any information about the origin of the star. 

Grenon l|1987l> proposed that, even if the present stel- 
lar position cannot be taken as an indicative of its birth- 
place, the average radius of its orbit, i? m , is kept close to 
the initial galactocentric radius of the stellar birthplace. 
This would make possible the use of R m in the derivation 
of radial constraints, like the abundance gradients and ra- 
dial variation of the age-metallicity relation. 

Wielen, Fuchs & Dettbarn 1)1996(1 criticized this con- 
clusion, pointing that R m must also be considerably mod- 
ified due to the orbital diffusion. Due to irregularities pro- 
duced by, for instance, giant molecular clouds and spi- 
ral arms, the average galactocentric radius of a stellar 
orbit R m , at time i, would not be directly linked with 
the galactocentric radius where the star was born, i? m (0). 
Wielen et al. have described the star scattering process 
due to the orbital diffusion, and have shown that after 
some Gyr it is impossible to know where the star was 
born, from R m , with a accuracy better than some kpc. 
Also, to explain the high metallicity dispersion in Edv93 
data, these authors suggest that the Sun and some other 
nearby stars were born closer to the galactic center, and 
have migrated to their present location due to the orbital 
diffusion. Supposing a very well mixed ISM, the authors 
propose that the stellar birthplace, Ri, can be roughly 



estimated from the age and metallicity of a star, by the 
equation 

Ri -R Q ~ -ll[Fe/H] - 0.53i 9 + 0.6, (2) 

where Rq is the present solar galactocentric radius, and 
ig is the age of the star, in Gyr. According to Widen et 
al., there is a relation between Ri and R m , given by 

(Rm-R®) =OA(R i -R ). (3) 

More recently, Nordstrom et al. l|1999f) have argued 
that the predictions by Wielen et al. may be overestimated 
and that R m seems to be a good indicator of the stellar 
birthplace. 

Our sample can be used to test whether Ri and R rn 
present the same behaviour predicted by Wielen et al. 
(1996). In Fig. E|we show a plot with values for Ri and 
R m . The solid line indicates the predictions by Wielen et 
al. I|1996fl . The scattering is too large to allow a significant 
conclusion. Moreover, the small range of R m in our sample 
makes difficult the comparison with the theory. 

Nevertheless, if Ri (or R m ) are good indicators of the 
stellar birthplace, we must find a relation between these 
parameters and the stellar age, since a finite time is nec- 
essary for a star to be able to move from its birthplace 
to the present galactocentric radius of the Sun. Thus, we 
expect to find stars born at several galactocentric radii 
only amongst the oldest stars, since the youngest would 
not have had time to migrate radially 1 or 2 kpc from its 
initial position. 

In Fig. we show how these parameters are linked 
with age. Panel a presents R m as a function of the stel- 
lar age, while panel b presents Ri. The data are for our 
sample. R m presents a well-marked correlation with age. 
It is clear that young objects have all R m « Rq, and that 
older stars present a higher proportion of objects coming 
from different galactocentric radii. If we assume that Ri 
represents the stellar birthplace radius, we do not find the 
same behaviour. It is difficult to understand how a star 
with less than 1 Gyr could come from galactocentric radii 
much smaller than Rq. Panels c and d present the same 
comparisons, for the sample by Edv93. The same conclu- 
sion can be taken from them. The comparison for Ri is 
even worse than for our sample. 

This figure shows that Eq. © is not suitable to cal- 
culate the stellar birthplace radius. However, it does not 
imply that the outwards migration of the Sun has not oc- 
curred. The formalism by Wielen et al. is developed for a 
well-mixed interstellar medium, in which the metallicity 
dispersion seen amongst the stars is produced by the or- 
bital diffusion. However, if there is instead a non-negligible 
metallicity dispersion in the gas, at each galactocentric ra- 
dius, the straightforward use of the stellar age and [Fc/H] 
in Eq. will only give a rough estimate for the birthplace 
radius. 

We have also tested this idea with a simulation of the 
stellar orbits diffusion. We have inserted tracer stars in 
an exponential disk under the 3D model for the Galaxy, 
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Fig. 7. Correlation between two proposed birthplace indicators. Panel a, R m versus age for our sample; panel b, Ri 
versus age for our sample; panel c, R m versus age for the sample of Edv93; panel d, Ri versus age for the sample of 
Edv93. The dashed line indicates the present galactocentric radius of the Sun. 



as described in Flynn et al. I|1996|l . New stars are born 
each time step, simulating a constant SFR. The stars are 
given small random kicks in their velocities at each time 
step. Stars near the solar circle get sufficient kicks for be- 
ing heated up to have values of the velocity dispersion 
ellipsoid similiar to the observed one. Closer to the galac- 
tic center, the heating is greater, increasing exponentially 
inward with the same scale length as the disk light distri- 
bution (circa 4 kpc), with less heating in the outer disk. 
Inside 2 kpc, the heating is cut off, to simulate less molecu- 
lar clouds there and to simplify the computation. Indeed, 
such stars are not expected to be deflected onto orbits 
which eventually take them out to the solar circle over 
the disk lifetime. 

At the end of the simulation, we find all stars which are 
on or near to the solar circle. We have computed the mean 
guiding radius of their present orbit, R m and compared it 
with the a priori known birth radius Ri. The relationship 



found between them is ((Ri — Rq)) ~ 2((R m — Rq)) as 
expected by the theoretical considerations performed by 
Wielen et al. l(UM|l . 

The problem does indeed come when one tries to re- 
cover the relationship from the data. The only way we 
have of recovering Ri from the data is to use [Fe/H]. The 
disk gradient is only about 0.1 dex/kpc, and the scatter in 
the local metallicities is (intrinsically) of the same order 
or larger (e.g., Edv93; Paper I). Here, by scatter, we mean 
the width of the [Fe/H] distribution. In our sample, i? m 
lies in the range 6 to 10 kpc, so that the birth radii of the 
stars lie in the range 7 to 9 kpc. Over such a short range in 
the disk, the mean mctallicity of the stars changes by only 
0.1 dex in either direction, whereas the intrinsic dispersion 
(over all ages) is 0.12-0.3 dex. 

Our results show that R m can surely give some infor- 
mation about the birthplace radii of the stars. But this 
information can be used only in a statistical sense, ac- 
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Fig. 8. Comparison of constraints for two samples, one internal to the solar radius, and the other external to it. Panel 
a, age-metallicity relation; panel b, metallicity distribution; panel c, age distribution. 



cording to Eq. ©, due to the intrinsic scatter in the R m 
relation introduced by the orbital diffusion. This must be 
kept in mind, when interpreting the results from the next 
section. 

6. Radial constraints to chemical evolution models 

Edv93 have used R m to investigate the radial variation of 
some chemical constraints, like the abundance ratios and 
the age-metallicity relation. We decided to use our data 
sample to derive constraints to the age-metallicity rela- 
tion, the radial variation of the metallicity distribution, 
the star formation rate and the abundance gradient in the 
disk. 

Due to the narrow range of R m in our data, we have 
divided the sample in only two bins according to their 
galactocentric radius: stars with R m < 8 kpc and stars 
with R m > 8 kpc, which corresponds to stars born in- 
ternally and externally to the solar galactocentric radius, 
respectively. These groups are composed by 158 and 222 
stars, respectively, adding up to 380 objects. 



Ideally, it should be possible to derive such constraints 
from a large data sample as ours. However, a close look 
at Fig. [7Ji shows that there is an intrinsic bias in the R m 
distribution. The younger objects were born in the vicin- 
ity of J?©, and are badly represented for R m ^ Rq. Any 
attempt to derive radial constraints are affected by this 
bias. 

Let us consider three possible constraints: the age- 
metallicity relation, the radial variation of the metallic- 
ity distribution, and the star formation history (from the 
age distribution), shown in Fig. [H] At first glance, there 
is no radial variation in the age-metallicity relation (Fig. 
|SJi), in opposition to the conclusions by Edv93, that the 
most inner parts of the Galaxy have evolved more rapidly. 
However, from Fig. |Ht we see that the sample internal 
to the solar radius does not comprise many young stars, 
since that have had not time to migrate to the solar vicin- 
ity. Since the younger stars are generally richer, the ra- 
dial variation presented in this figure cannot be consid- 
ered real. This explains why the metallicity distribution 
of the inner sample presents a paucity of metal rich stars, 
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compared to the metallicity distribution of the outer sam- 
ple (see Fig.|SjD). The different age distributions in Fig. [HJ; 
also must not be interpreted as radial variation of the star 
formation history, because our sample does not have an 
homogeneous R m distribution as a function of the stellar 
age. 

In Fig. we show the metallicity as a function of the 
birthplace (upper panel). The different symbols separate 
the objects according to their age, namely: stars with ages 
lower than 1 Gyr (filled squares), stars with ages of 1-3 
Gyr (open circles), stars with ages from 3 to 6 Gyr (dotted 
triangles), and stars with ages greater than 6 Gyr (upside 
down triangles). We can see that the metal-rich objects are 
those that were born in the vicinity of the Sun, whereas 
the metal-poor objects come from several galactocentric 
radii. The same behaviour is present in Edv93 data (lower 
panel). 

The present sample is not particularly useful in the 
study of radial abundance gradients, in view of the un- 
certainties in the determination of i? m and the magnitude 
of the gradients, which is comparable to the average dis- 
persion in the abundances as taken in the whole range of 
galactocentric radii. However, it is interesting to note that 
an average [Fe/H] gradient of -0.07 to -0.09 dex/kpc as 
adopted by Maciel ( 2002) for open cluster stars fits rather 
nicely our data points with ages under 6 Gyr as shown by 
the filled squares, open circles and dotted triangles in the 
upper panel of Fig. |5J 

The conclusion by Edv93, that the inner parts of the 
Galaxy have experienced a more rapid enrichment, seems 
somewhat premature. There is no unbiased coverage in 
R m and [Fe/H] in their sample to allow such conclusion, 
just like in our sample. Although radial constraints are 
extremely needed for the understanding of the chemical 
evolution of the Galaxy, presently there is no sample of 
late-type dwarfs that can be used in their derivation. 

7. Implications for Theoretical Models 

Our results for the AVR, while fairly consistent with in- 
dependent work, add to the growing evidence that heat- 
ing by scattering by GMCs, spiral structures, and accre- 
tion, are insufficient to explain the observations. Here we 
summarize the problem, and suggest that a different in- 
terpretation of the AVR should be re-investigated. In the 
following, we mostly consider velocity dispersions with re- 
spect to a cylindrical coordinate system, where r, 6 and z 
represent the radial, azimuthal and vertical components of 
the space velocity, respectively. For Solar neighbourhood 
stars, we can adopt a w ja u = a z ja r . 
Our major results are: 

1. Excellent power law behavior of all the velocity dis- 
persion components with age, with an exponent of 
about x = 0.24-0.35, similar but somewhat smaller 
than Wielen's ^Wj\ original estimate and some more 
recent results (Jahreifi & Wielen 1997] Fuchs et al. 
I2UUT)) based on HIPPARCOS data. Our exponents are 
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Fig. 9. Metallicity-i? m diagram. Upper panel: our data; 
lower panel, Edv93's data. The symbols indicate: filled 
squares, stars having ages lower than 1 Gyr; open circles, 
stars aging 1-3 Gyr; dotted triangles, stars having ages 
from 3 to 6 Gyr; upside down triangles, stars with ages 
greater than 6 Gyr. 

increased by about 0.1 if \W\ weighting is used, but 
there is much more scatter in the fits for each com- 
ponent as well as other problems described above. We 
note that Jahreifi & Wielen <|1997[) and Fuchs et al. 
(2001) used \W\ weighting. 

2. We also find that a w reaches values around 30 km/s 
for the oldest stars, while the total velocity dispersion 
reaches 70 km/s. 

3. Another constraint on theories for disk heating is the 
rate of axial ratios of the velocity ellipsoid, for which 
we find o~ w : o~ u ~ 0.45 to 0.55 in most age bins, with a 
spread of about 0.1 between the different age groups, 
and little dependence on age except for the two small- 
est age groups. These values are a little smaller than 
previously adopted values of 0.5-0.6 (e.g. Wielen lT977|l . 
but the lack of age dependence agrees with previous 
work (Dehnen & Binney 1998). 

These results appear inconsistent with calculations 
based on scattering from GMCs and spiral waves, pri- 
marily because none of these calculations can give a large 
enough velocity dispersion for old stars. The other con- 
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straints are more difficult to assess. Lacey l|1984|l extended 
Spitzer & Schwarzschild's ifllMl IT§5l|) Fokker-Planck for- 
mulation of the problem to 3-dimensional disks and found 
x w 0.25, a z ja r ~ 0.8, suggesting that the velocity el- 
lipsoid was too round. This led to a revival of the idea 
that transient spiral arms could heat the disk. Through 
resonances of stellar orbits with wave frequencies, primar- 
ily radial heating is produced, with some of it channeled 
into the vertical direction. Sellwood & Carlberg ©84), 
Carlberg & Sellwood <fT5%5)l . and Carlberg fTSE^ all found 
that transient spiral arms can give <J z /a r sa 0.5, but ob- 
tain inefficient vertical heating, with x ~ 0.2. This led 
Jenkins & Binney (19Q and Jenkins l|1992fl to analyze 
the problem by combining scattering from clouds and spi- 
ral waves. They found that clouds are indeed more effi- 
cient vertical heaters in the presence of spiral waves, but 
still obtained x < 0.3. Jenkins {TO92) proposed that the 
addition of disk accretion to the scattering mechanisms 
could alleviate the problem, giving x « 0.4, but obtained 
a z /cr r in 0.65, which appears too round (especially com- 
pared to our results). Furthermore, it is not possible to 
reach vertical velocity dispersion of ~ 30 km/s if young 
stars have a w = 2-5 km/s (Asiain et al. 119991 see Kroupa 
EZUOH Fig. 7). Kroupa l|2"U0"2"jl investigated the addition of a 
component due to unbound star clusters which expand at 
virial or super-virial velocities following fast gas removal. 
Such a component adds interesting effects to the tail of 
the stellar velocity distribution, but the cr^-age relation 
requires a star formation rate that has been strongly de- 
clining with time or a cluster mass function which has been 
shifting to less massive clusters with time to the present 
epoch. 

Perhaps the discrepancy is due to the many approxi- 
mations made in the Fokker-Planck theory, which are well- 
summarized by Villumsen ( 1985). Villumsen presented nu- 
merical simulations that overcame many of these short- 
comings. After a brief initial period of fairly rapid in- 
crease of heating, the velocity dispersions flatten out to 
= 0.25 (in agreement with Lacey's analyt- 
ical treatment), and x w « 0.31. By the end of the sim- 
ulations (~ 5 x 10 9 yr) the ratios of velocity dispersions 
were ay : ag : a x = 1.0 : 0.7 : 0.6. Thus the velocity dis- 
persion is too round, although not as round as found by 
Lacey's (|1984fl analytic result. Judging from the papers 
cited above, transient spiral arms would reduce the a z /a r 
ratio, but the addition of disk accretion should increase it 
again, by an unknown amount. Extrapolating Villumsen's 
velocity dispersions to 10 10 yr indicates that the total dis- 
persion will not exceed 40 km/s, while a z will only reach ~ 
15 km/s, compared to about 70 and 30 km/s, respectively, 
found here for the oldest stellar group. A major approx- 
imation made by Villumsen ljl985[) is that the number 
of clouds in the simulations was greatly reduced to 200 
for computational reasons while the masses of the clouds 
were reduced by the square root of this factor to compen- 
sate. It is not known how the diffusion argument on which 
this scaling is based might affect the results. However, re- 
cent high-accuracy simulations of 1O 6 M GMC heating by 



some of us (Hanninen & Flynn 20Q]2|) give an even weaker 
increase of the velocity dispersion with time, although the 
velocity ellipsoid is flatter. 

The situation may actually be worse than this. In a 
pair of little-cited papers, Yasutomi & Fujimoto (1989 
1991) presented a significantly more sophisticated treat- 
ment of the star-cloud scattering problem, including the 
non-Keplerian potential of individual clouds (assumed 
Plummer), the nonlinear part of the background disk 
gravitational field (neglected in earlier work which as- 
sumed harmonic oscillation through the linear epicyclic 
approximation), the galactic shear of stars and gas clouds, 
and the non-uniform distribution of massive molecular 
clouds in the Galactic plane (annulus-like 4-8 kpc ring and 
spiral-like distributions), but restricted to two dimensions. 
The main approximation made is that stars only interact 
with their nearest gas cloud. The major effect discovered 
was that statistically anisotropic encounters (direct and 
catching-up encounters) lead to important systematic ac- 
celerations and decelerations of the stars, corresponding to 
dynamical friction. This leads to very slow growth of the 
velocity dispersion for times > 1 Gyr, with x m 0.20 for a 
uniform distribution of clouds, and even smaller values ~ 
0.10 to 0.15 for clouds distributed in an annulus or spiral 
arms. In all cases the maximum total velocity dispersion 
was only 30 to 50 km/s, a value reached quickly (< 1 Gyr), 
before dynamical friction effects become important. These 
results led Yasutomi & Fujimoto to suggest that if a scat- 
tering process is to account for the observations, either 
gas clouds must have been more massive and numerous 
billions of years ago fSemenzato ll987)l . and/or there must 
exist a high surface density of massive dark objects, such 
as the 10 6 M Q black holes proposed by Lacey & Ostriker 
( 1985J) . The generalization of these calculations to three 
dimensions by Yasutomi & Fujimoto Ijl991| showed that 
the general results are robust. In none of their models did 
the total velocity dispersion reach values in excess of 35 
km/s, and in fact slowly decreased with time due to the 
dynamical friction effect, after the initial transient heating 
episode for models with a uniform spatial distribution of 
clouds. Furthermore, they quote a result for the velocity 
ellipsoid oy : ag : a z = 1.0 : 0.7 : 0.7, clearly too round 
compared to the results found here or in earlier work. 

Kokubo & Ida 1 1992(1 reexamined the problem using 
direct integrations of epicycle orbits to obtain probabil- 
ity distributions for the change in orbital parameters, but 
obtained results similar to the other methods. After an 
initial rapid increase, the later phases gave x = 0.25 for 
all three velocity dispersion components (smaller than our 
result) and a z /cr r ss 0.6 (slightly larger than our result). 
For GMC masses of 10 6 M Q , the total velocity dispersion 
only reaches about 30 km/s while a w only reaches 10 km/s 
after 10 10 yr. Kokubo & Ida find that if the GMC mass 
is increased to 1O 7 M , cr tot increases to about 50 km/s, 
while a z reaches about 20 km/s. Kokubo & Ida's claim 
to the contrary, these values are still significantly smaller 
than found here and elsewhere. 
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Thus the current situation is that none of the above in- 
vestigations of scattering/heating models can account for 
either x, a z j a r or the total or vertical velocity dispersions 
found here and elsewhere for older stars, even when shear 
motions are included. Including more physical realism has 
apparently only exacerbated or accentuated the problems 
in earlier work. 

The ratio a z ja r has also been estimated for at least 
two other disk galaxies, NGC 488 (Gerssen, Kuijken & 
Merrifield ITWHl and NGC 2985 (Gerssen, Kuijken & 
MerrineldEUnni)- The values that they found (0.70 ± 0.19 
and 0.85 ±0.1, respectively) are larger than found in the 
solar annulus, although still significantly less than unity. 
The existence of a broad range in this ratio among disk 
galaxies would significantly complicate the interpretation, 
especially since there are several different mechanisms for 
producing a ratio of a given value and uncertainty in the 
result for each mechanism. We therefore concentrate on 
the absolute values of the velocity dispersions themselves, 
and confine our interpretation to the solar annulus where 
the data is probably most accurate. 

One possibility is that the lifetimes of the clouds were 
assumed to be too large (infinite). Fujimoto (1980}, using 
a significantly different stochastic approach to the problem 
(theory of random walks in the field of a harmonic oscil- 
lator), assumed that GMC lifetimes are small, ~ 10 7 yr. 
At the time of this work, the large lifetimes (~ 3 x 10 s yr) 
suggested by Solomon, Sanders & Scoville (|1979[) . based 
on coagulation requirements, were popular, and so sev- 
eral papers dismiss Fujimoto's calculations as unrealis- 
tic because of the small adopted cloud lifetimes. However 
subsequent work has strongly suggested smaller lifetimes 
about equal to what Fujimoto adopted for "CMC's (see 
Elmegreen 2000 Hartmann, Ballesteros-Paredes & Bcrgin 
(200IJ, and references therein). The effect of the short 
lifetimes is significant because it provides a strong fluc- 
tuating force component which is missing in all other 
work. Kroupa (2002), unaware of Fujimoto's work, inde- 
pendently pointed out how the small cloud lifetimes would 
add a fluctuating acceleration to stars because they expe- 
rience a time-varying potential on the timescale of their 
passage by each cloud. However, Fujimoto obtained to- 
tal velocity dispersions and vertical velocity dispersions of 
only 40 km/s and 6 km/s after 10 10 yr. This result must be 
partly related to Fujimoto's small assumed GMC mass of 
2 x 10 5 M Q , but the scaling of the heating rate with GMC 
mass derived by Kokubo & Ida ill 992 i suggests that even 
this short-lifetime effect will be insufficient to provide the 
large velocity dispersions of the older stars. 

In order to use heating to account for the AVR, it may 
be necessary to invoke very massive black holes (Lacey & 
Ostriker 1985). Simulations by Hanninen & Flynn (2002) 
show that if the total mass of the halo consists of 1O 7 M 
black holes, the disk could be heated from a present veloc- 
ity dispersion of 18 km/s to 80 km/s in about 8 Gyr (i.e., 
for an 8 Gyr old population), but if the present velocity 
dispersion is 10 km/s, the disk would only be heated to 
about 60 km/s even for 15 Gyr old populations. Only a 



combination of GMCs and 10 7 Mq black holes comprising 
half the halo could heat the disk sufficiently in 12 Gyr, but 
in that case a z /a r w 0.63, somewhat larger than found 
here, although marginally consistent with other results. 

Another mechanism which could plausibly result in a 
velocity dispersion-age relation is the stochastic heating 
due to minor mergers of the Galactic disk with satellite 
galaxies or passage of dark matter substructure through 
the disk. The dynamical friction force exerted by the disk 
converts a fraction of the satellite kinetic energy into ver- 
tical disk energy. Older stars should have suffered more 
events, either from more satellites and dark matter sub- 
structure, or repeated orbital passages of an individual 
object until it is tidally stripped, resulting in a velocity 
dispersion-age correlation. This process has been stud- 
ied analytically and numerically by a large number of 
authors (see Velazquez & White HMfl Font et al. I2UUT1 
Taylor & Babul 2001J Benson et al. 120031 and references 
therein). Unfortunately most of this work was aimed at 
the problems of disk disruption, and then matching the 
scale heights of present-day disk galaxies. The problem 
is made more difficult by the uncertain partitioning of en- 
ergy deposited in the disk between disk heating and global 
modes such as warping and bars (which themselves could 
heat the disk). Font et al. ( 2001) used simulations of disks 
heated by the ensemble of subhalos in CDM model halos 
to infer that the disk heating rate (as measured by an es- 
timated diffusion coefficient) is less than the rate required 
to match observations of the solar neighborhood, although 
the discrepancy is not very large, and depends somewhat 
on the adopted distribution of satellite orbits and espe- 
cially the mass distribution of the CDM subhalos. Benson 
et al. H20p3) used similar ensembles of halos with a wide 
range of masses to show that the predicted scale height 
distribution of present-day models (i.e. after billions of 
years of heating) compares favorably with the distribu- 
tion of scale heights determined by Bizyaev & Mitronova 
( 2002) for a complete sample of 60 edge-on galaxies, but 
only when some molecular cloud scattering is included to 
avoid a low-scale height tail of the distribution. However 
they did not give information concerning how the velocity 
dispersion increases with time, its absolute value, or the 
ratio of vertical to radial velocity dispersion. After 4 Gyr 
of the orbit of a single satellite, the disk vertical energy in- 
creases by about 10% to 80% for various models of Benson 
et al. H2003[) . As a velocity dispersion increase (i.e. square 
root of energy increase), this is very small compared to 
the factor of 2 or so increase that the Milky Way exhibits 
over a similar time, but the number of discrete satellite 
bodies could be large, at least initially (about 100 in the 
Font et al. models). 

One problem we see with this model is that tidal strip- 
ping by the disk continually decreases the flux of bom- 
barding subhalo structures, so one expects that eventually 
the velocity dispersion-age relation should flatten out at 
ages less than some value, or at least show a change of 
slope due to the decreasing effectiveness of this process 
at more recent times. Such a flattening is seen in the ve- 
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locity dispersion-age relations shown in Font et al. I|2001[l . 
No such feature is seen in the observations presented here, 
supporting their conclusion that heating by merger with 
subhalos is unlikely to account for the relations. 

We cannot rule out the possibility that minor mergers 
could have imprinted sudden changes in the slope of the 
AVR, especially given the cumulative observational evi- 
dence that the Milky Way faced and is presently facing 
mergers with some dwarf satellites (e.g., Ibata, Gilmore & 
IrwinHnHl Gilmore, Wyse & Norris EStMl Newberg et al. 
ISUTEfl Yanny et al. l2HIT31 Rocha-Pinto et al. l2UU31 Majewski 
et al. l2003[l . It is tempting to ask whether slope changes in 
the AVR could be used to recover the past merging history 
of the Galaxy. Mergers could also trigger star formation 
bursts in the disk, as Paper II suggests. Nevertheless, we 
have found no clear correlation between the star formation 
history (derived in Paper II) and the AVR. 

Similarly, it appears that vertical heating of a galactic 
disk by bending instability (Sotnikova & Rodionov 2003), 
while possibly important due to the relatively long sat- 
uration timescale of the instability, appears to give sig- 
nificant "shelves" in the velocity dispersion-time relation. 
These features appear at the time of maximum ampli- 
tude of the nonaxisymmetric and then symmetric bend- 
ing, at relatively small times (0.5 to 1.5 Gyr in the mod- 
els of Sotnikova & Rodionov), followed by a flat velocity 
dispersion-time relation. Again, such saturation is not seen 
in the observations, although it is uncertain if the mecha- 
nism, in conjunction with the other mechanisms discussed, 
could avoid the saturation. However the radial velocity 
dispersion component actually decreases in time due to 
this mechanism, again suggesting that it may be impor- 
tant, but would need to be combined with other mecha- 
nisms to explain the observations. The results for heat- 
ing by the bending of a bar presented by Sotnikova & 
Rodionov are similar (except that the radial dispersion in- 
creases somewhat with time) , with a factor of 2-3 increase 
in the vertical dispersion produced in about 1.5 Gyr. 

A very different, but entirely plausible, way for dissi- 
pative gas systems to account for the AVR that was first 
proposed by Tinsley & Larson l|1978[) . They showed, in 
the context of Larson's (1976) two-fluid models for disk 
galaxy formation and evolution, that the gradual dissipa- 
tion and settling of the gas into a thinner and cooler layer 
naturally results in a reduction of the vertical velocity dis- 
persion from about 50 km/s at age 3 Gyr when the disk 
becomes well-defined, to about 15 to 12 km/s for ages be- 
tween 10 and 15 Gyr. The stellar velocity dispersion just 
reflects the evolution of the combined internal velocity dis- 
persion and (for the older stars) the vertical motion of the 
gas disk. 

Advantages of the Tinsley & Larson dissipation model 

are: 

1. The sign of the a(t) correlation and perhaps the ve- 
locity dispersion for old stars are accounted for fairly 
naturally, without any fine-tuning of parameters (for 



later work, see below) or introduction of, for example, 
extremely massive scattering entities; 
2. The correlation between velocity dispersion and aver- 
age metallicity, both of which should depend on age, 
can be accounted for essentially using only one phys- 
ical process, since turbulent dissipation and the star 
formation rate are closely coupled in disk evolution 
models. Disk heating mechanisms that involve stellar 
scattering are decoupled from chemical evolution and 
the star formation rate. 

The problem with the model is that it is not clear what 
the predicted functional form of the AVR and the corre- 
sponding components should be. In addition, the evolu- 
tion of the gas velocity dispersion is controlled in part by 
the adopted parameterization of the star formation rate in 
terms of density or other variables, which is problematic 
because of the notorious uncertainty in such parameteri- 
zations. 

Later fluid models for the settling of initial hot proto- 
galactic gas into the equatorial plane have become more 
detailed in input physics making use of the chemodynam- 
ical prescription of galaxy evolution, including different 
forms of star-formation feedback, a number of heating 
and cooling processes, and an accounting of the multi- 
phase character of the interstellar medium. Although more 
parameters are naturally required to accommodate the 
increased number of processes, their choice is not arbi- 
trary but formulated from plasmaphysical and astrophys- 
ical studies. It can also be demonstrated that even if the 
parametrization of the star-formation rate is taken arbi- 
trarily in the chemodynamical prescription it will install a 
dependence on density and temperature of the warm gas 
phase in a very narrow range because of self-regulation 
processes (Koppen et al. 19951 IT9"H%j) . Burkert, Truran & 
Hensler (1992) provide a useful example, although only the 
vertical dependence is included in their one-dimensional 
fluid model. The results are very similar to Tinsley & 
Larson, with the velocity dispersion decreasing from ~ 
40 km/s to ~ 20 km/s after about 5 Gyr. Although they 
ascribe this decrease to a decline in the star formation 
rate (hereafter, SFR) by a factor of 10 due to gas de- 
pletion (they assumed the SFR varies with the square of 
the gas density), this is equivalent to Tinsley & Larson's 
gaseous dissipation effect, since it is the SFR that drives 
the "cloud" or "turbulent" motions of the gas from which 
the stars form in both cases. Besides the advantages men- 
tioned above, the work of Burkert et al. ( 1993J) shows that 
their models can account for the observed vertical distri- 
bution of stellar density and velocity dispersion. 

An extension of this chemodynamical model to two di- 
mensions and more physical processes has been presented 
by Samland, Hensler & Theis l|1997j> . A notable result of 
the Samland et al. models (their Fig. 3) is that, although 
the total (halo + bulge + disk) SFR has an early peak at 
~ 2 Gyr and then declines by a factor of five or so, the 
disk SFR is actually fairly flat. As explained by Samland 
et al., this is because the disk has only become a relatively 
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redshift properties, and have used SPH + N-body meth- 
ods rather than fluid models, making it difficult to esti- 
mate velocity dispersions at late times when small SFRs 
mean that the number of newly-formed star particles is 
small. However, the velocity dispersion behavior of such 
models, when available, are very similar to those discussed 
above. For example, Steinmetz & Miiller (1095) state that 



Fig. 10. Temporal evolution of the vertical velocity dis- 
persion of the cloudy gas phase in the solar vicinity (at 8.5 
kpc galactocentric distance) as the result of the chemo- 
dynamical Milky Way model by Samland et al. ©97), 
adapted from Samland l|1995[) . 



discrete entity after about 7-8 Gyr, while the total dura- 
tion of the model is 15 Gyr. This slow disk formation in 
their models can be traced back to the large cloud collision 
(i.e. dissipation) timescales adopted. Even beginning at 7 
Gyr, the disk SFR appears to only decrease by a factor 
of three or so, suggesting that the velocity dispersion will 
not decrease much with time. 

Fig. 1101 shows the time history of the vertical velocity 
dispersion of the cloudy gas phase in the solar vicinity (at 
8.5 kpc galactocentric distance) obtained in the chemo- 
dynamical Milky Way model by Samland et al. 1)1997(1 . 
Since stars are formed from the cool phase this also rep- 
resents the velocity dispersion of the newly formed stars. 
One can discern that the vertical dispersion has early os- 
cillations between 45 and 65 km/s, when gas starts to 
accumulate in the equatorial plane and partly cools dissi- 
patively but is still heated due to the continuing collapse. 
Then, after the disk has formed about 8 Gyr ago (in the 
model), the vertical velocity dispersion declines, flatten- 
ing out at about 25-30 km/s over the past 5 Gyr. For the 
comparison of this value with observed solar vicinity stars, 
one has to realize that the result is derived as an average 
over a few vertical grid cells and thus only represents a 
thick disk component of almost 700 pc height because of 
its coarse spatial resolution. The ratio a z ja r is about 0.5 
at most radii in the late-time model, in good agreement 
with our empirical result. 

Other models for disk galaxies that include dark mat- 
ter and cosmologically-motivated initial conditions have 
tended to concentrate on chemical evolution and/or high- 



for stars in age intervals [3,8.5], [1,3], and (indirectly) < 1 
Gyr, the velocity dispersions fall from 80 to 40 to 20 km/s, 
respectively, at R ~ 10 kpc. Their Fig. 8b suggests that 
at late times oy > <jq a z over most of the model galaxy 
that is relevant to the solar neighborhood, but small num- 
ber fluctuations and artificial heating by scattering from 
the model dark matter particles make the results difficult 
to interpret. Once again, the agreement with the density 
distribution and chemical properties (e.g. metallicity gra- 
dients, see Steinmetz & Miiller 1994) of real galaxies makes 
this cooling interpretation of the a(t) correlation attrac- 
tive from the point of view of "Occam's razor." 

More recent SPH/N-body models of disk galaxies by 
Berczik ( 19991 l2~()()()|l . which include more stellar compo- 
nents and different parameterization of stellar energy feed- 
back, as well as a different star formation rate prescription 
(basically a virial theorem condition on the SPH parti- 
cles rather than either a density dependence or a density 
threshold) do not give the time-dependence of the stel- 
lar velocity dispersion, but indicates that final values are 
a z w 20 km/s over the whole disk, with <j r > a z near the 
sun (see Berczik Figs. 6-8). 

Thus the time dependence of the disk velocity depen- 
dence and its absolute scale seems robust over a fairly 
wide range of independent models, agreeing overall with 
Larson's l|1976f) early work. It is unfortunate that the rel- 
evant kinematic information from the detailed disk evo- 
lution models of Raiteri, Villata & Navarro l|1996p and 
Fricdli & Benz (1995) are not available to test this robust- 
ness further. The published star formation rates suggest 
that the behavior will be similar, although there are larger 
apparent temporal fluctuations in the late disk SFRs and 
presumably in the stellar velocity dispersions. 

However unlike the scattering models, which have 
great trouble accounting for the large velocity dispersion 
of older stars, the cooling models have a problem in pro- 
ducing small enough velocity dispersions at late times: 
Instead of power law cooling, the velocity dispersion ap- 
proaches a nearly constant value. Larson's early models 
apparently gave the smallest present day velocity disper- 
sions of 12-15 km/s, while all the other models discussed 
above tend to flatten out at about 20-35 km/s. In fact 
Tinsley & Larson l|1978f) and Larson Q1979[) were aware 
of this problem, and suggested that turbulent dissipation 
within the disk layer (rather than settling of the disk itself 
on larger scales) could lead to the requisite small veloc- 
ity dispersion. The reason this does not occur in the disk 
formation models is that at late times the SFR is small, 
so that the star formation uses available gas on a long 
timescale, leading to a relatively constant SFR. Since the 
models tend toward self-regulating equilibrium, with dissi- 
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pation balancing SF (see Koppen, Theis, & Hcnslcr 1995 
for one-zone models), the gas and hence stellar velocity 
dispersion cannot decrease. 

Nevertheless, one reason for the apparent large veloc- 
ity dispersions in all the disk formation models at late 
times might be lack of spatial resolution, such that the 
calculation is yielding an average velocity dispersion over 
a height much larger than that of the thin disk. 

A likely more important aspect of these models is that 
they all use some form of sub-grid (unresolved) dissipation 
whose timescale is essentially a parameter in the models. 
For example, in the chemo-dynamical collapse models, the 
dissipation comes through cloud collisions, but the dissi- 
pation timescale is set to about an order of magnitude 
larger than the cloud collision timescale based on the ge- 
ometric cross section. This allows for slow early collapse, 
but also forces the gas to resist a reduction in velocity dis- 
persion below about 20 km/s. It is relatively easy to show 
using one-zone analytic models that the self-regulated gas 
velocity dispersion will scale with the SFR as (SFR) 1 / 3 
- see Scalo & Chappell <|1999[) . In fact a number of in- 
dependent simulations of interstellar MHD turbulence in- 
dicate that the dissipation time is rapid, of order a cross- 
ing time (Mac Low et al. 119981 Mac Low 119991 Ostriker, 
Gammie & Stone [1999 ), even in the presence of stellar 
heating sources (Avila-Reese & Vazquez-Semadeni^OOl). 
Thus we ascribe the inability of collapse models to settle 
to a small-enough velocity dispersion to an underestima- 
tion of the dissipation rate. If the dispersion were really 
set by an equilibrium between SF and dissipation, this 
model would be unable to account for the continued de- 
crease in velocity dispersion at relatively recent times (< 
1 Gyr). We think the problem is not serious, since, mov- 
ing away from one-zone models, more detailed turbulence 
simulations show that the turbulence dissipates rapidly in 
spatial regions between star-forming events (Avila-Rccsc 
& Vazquez-Semadeni .2001}, and new stars will form from 
this gas. 

In fact it is obvious that this model has no problem in 
accounting for the low velocity dispersion of young stars, 
since the velocity dispersion of the interstellar gas is ob- 
served to have about the same value as that for the stars. 
The only problem involves understanding why the exist- 
ing simulations do not dissipate turbulence sufficiently to 
match the observed present-day gas velocity dispersions. 

In the future it may be possible to empirically dis- 
tinguish between the dissipation and heating models by 
comparing the models with observations of a z /a r in disk 
galaxies. Shapiro, Gerssen & van der Marel <|2003[) have 
used major- and minor-axis kinematics to estimate a z /a r 
for a number of disk galaxies of various Hubble types. 
The resulting values of er 2 /oy are in the range 0.5 to 0.8 
with considerable uncertainties, and are consistent with 
no trend with Hubble type at the lex level, or with a 
marginally significant trend of decreasing a z /a r with later 
Hubble type. 

Unfortunately, it is not yet possible to test these re- 
sults with predictions of the two models. Simulations by 



Jenkins & Binney l|1990fl seem to predict a decrease with 
increasing molecular gas fraction and spiral arm strength, 
so with advancing Hubble type. The uncertainties in the 
empirical ratios are large enough that such a trend cannot 
be uncovered at the la level. The galaxy collapse models 
can produce cr z /a r ratios in the observed range, but the in- 
ability of the models to predict the ratio at late times, and 
the lack of modeling of disks of different Hubble types, pre- 
clude a prediction of any trend. Hopefully further models 
of disk kinematics in galaxy formation models that result 
in a range of Hubble types will yield such information. 
One trend that could in principle be tested is that a z /a r 
varies from 0.6 in the dense inner regions to 0.48 in the 
outermost region in the model of Samland et al. I|1997[) . 
It is not known whether this variation is generic for the 
dissipation model. 

8. Conclusions 

Following our former investigations on the age-metallicity 
relation and the star formation history, we have used the 
chromosphcric ages for a sample of late-type dwarfs to 
investigate the chemokinematical properties of the solar 
vicinity. The use of chromosphcric ages warrants an al- 
most complete coverage of the galactic disk history with a 
sample consisting of stars in the same colour range, some- 
thing that cannot be achieved from other methods. 

The exponents for the AVR that we have found frus- 
trantingly seem to span the whole range of values which 
can be used to constrain the theory on the disk heat- 
ing. Nevertheless, after taking into account a statistically 
more meaningful binning and metallicity corrections to 
the chromosphcric age, we arrive at an exponent ranging 
from 0.26 to 0.31, for a tot oc (1 + t/r) x . 

We have shown that the vertex deviation of solar 
neighbourhood late-type stars varied irregularly between 
the present time and 10 Gyr ago. No clear dependence of 
the vertex deviation on the age of the stellar group was 
found in our analysis. The vertex deviation also varies sub- 
stantially when different selection criteria are applied to 
the sample. We give support to Palous ( 195rj|l idea that 
the vertex deviation is caused by the contamination of the 
sample by superclusters and moving groups stars. 

Our sample was used to test the idea that there is a sta- 
tistical relation between the average orbital galactocentric 
distance and the galactocentric radius where the star was 
born. While our simulations provide support to Wielen et 
al. (1996) results, we have found that local data does not 
have a good coverage in orbital galactocentric distance to 
allow a empirical testing of this relation. Nevertheless, the 
relation between R m and age that we have found supports 
Edv93's use of R m as a reasonable indicator of the stellar 
birthplace radius. 

Nevertheless, since young stars born at a galactocen- 
tric radius different from Rq have not had time to scatter 
into the solar neighbourhood at representative numbers, 
we believe that there is no possibility to find the radial 
variation of the age-metallicity relation or star formation 
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rate using local stars only. A much larger data sample, as 
that expected to be provided by GAIA, could solve this 
problem. 

A consideration of numerous results of theoretical 
studies of secular heating mechanisms, most involving 
gravitational scattering from large inhomogeneities in the 
gas or star system, indicates that none of these is capable 
of explaining satisfactorily the large velocity dispersions 
of the oldest disk stars. The problem is especially serious 
when one considers the dynamical friction effect found by 
Yasutomi & Fujimoto $T$M ES2) ■ 

Stochastic heating mechanisms involving satellite in- 
teractions are possible, but should (and do) give an AVR 
that flattens out for relatively young (few Gyr) stars, a 
feature that is not observed. Instead, we suggest that a 
very plausible alternative is the one suggested long ago by 
Tinsley & Larson l|1978fl . that the velocity dispersion of 
stars just reflects the velocity dispersion of the gas from 
which they formed, and that the turbulent motions of this 
gas has been declining continually since the formation of 
the disk. We cite a number of more recent disk evolution 
papers that find similar behaviour. 

While all the models successfully account for the dis- 
persion of the older stars, we are unable to say whether 
the simulations produce the low velocity dispersions ob- 
served for young stars, due to a combination of resolution 
effects and uncertainties in input physics. However this 
is obviously only a problem for the simulations, and not 
for the idea itself: the real observed interstellar medium 
does have a velocity dispersion that is similar to the veloc- 
ity dispersion of the youngest stars, so we are assured that 
the model must give the correct result at the present time. 
Hopefully future studies of disk evolution will improve the 
situation. 

There is a straightforward way to test and distinguish 
the heating and cooling models if stellar and gas velocity 
dispersions for face-on high redshift galaxies can be ob- 
tained in the future. The heating model predicts that the 
stellar velocity dispersion should be small at such early 
times, while the cooling models predict that the gas ve- 
locity dispersion should be large. We think it unlikely that 
both conditions can occur simultaneously. 
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